A Three-Dimensional Model of Small Signal Pree-Electron Lasers 
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Coherent electron cooling is an ultra-high-bandwidth form of stochastic cooling which utilizes the 
charge perturbation from Debye screening as a seed for a free-electron laser. The amplified and 
frequency-modulated signal that results from the free-electron laser process is then used to give an 
energy-dependent kick on the hadrons in a bunch. In this paper, we present a theoretical description 
of a high-gain free-electron laser with applications to a complete theoretical description of coherent 
electron cooling. 



I. INTRODUCTION 

Coherent electron cooling (CeC) [l| is a new cooling 
method for intense relativistic hadron beams, to be im- 
plemented first at the proposed MEcRHIC/cRHIC up- 
grade to the RHIC accelerator at Brookhaven National 
Lab. Schematically similar to the stochastic cooling al- 
ready implemented at RHIC 0], CeC has the advantage 
that its coherent bandwidth is on the order of the res- 
onance wavelength of the operating free-electron laser, 
so that the cross-correlation that leads to heating and 
therefore saturation of the stochastic cooling system is 
not encountered in CeC. 

To achieve a complete theoretical description of Co- 
herent Electron Cooling, models for the propagation of 
a phase space perturbation through thepick-up 0] and 
kicker [3| were developed and presented [g. All these cal- 
culations are based upon an infinite electron beam with 
K — 2 energy spread jl8|. However, an exact analytical 
solution for the high gain free-electron laser in the small 
signal regime, given an initial phase space perturbation, 
had not yet been developed. 

A number of analytical models have been developed for 
the transverse laser profile for an FEL. A set of equations 
for the full dynamics of a three-dimensional FEL with be- 
tatron oscillations were first written down in Q . Univer- 
sal scaling for the gain of the FEL in terms of the energy 
spread, emittance and focusing properties were developed 
in A fully three-dimensional Maxwell- Vlasov equa- 
tion was studied in Q and ultimately a procedure for ex- 
act and variational solutions to the laser eigenmodcs was 
presented in Q. These results focus primarily upon an 
eigenmode of the generated laser field, without consider- 
ation for developing solutions to the phase space density 
of the electron bunch. In this paper, we present a theo- 



* Electronic address; ! swebb@gr ad. physics. sunysb.edu] 
t Collider- Accelerator Department, Brookhaven National Labora- 
tory 

Electronic address: 



retical picture of the full dynamics of the electron phase 
space distribution, neglecting betatron oscillations, with 
an intent of using this result in application to CeC. 

In Section 2, we present an overview of the configu- 
ration of Coherent Electron Cooling, and discuss briefly 
the existing results in the pick-up and kicker sections. 
With the context of this work in mind, we then present 
a derivation for the dynamics of a high-gain free-electron 
laser seeded with an initial phase space perturbation in 
Section 3. This leads to an equation for an arbitrary 
transverse distribution of an otherwise infinitely long 
electron beam. In Section 4, we analyze the case of an 
infinitely wide beam, which leads to a Green function for 
the 3D FEL process with an infinitely wide beam. A 
mode expansion method is considered for a finite beam 
in Section 5. To conclude, we consider the specifics of 
applying these results to Coherent Electron Cooling in 
Section 6. 



II. OVERVIEW OF COHERENT ELECTRON 
COOLING 

Coherent Electron Cooling is schematically identical 
to stochastic cooling with a pick-up which gathers 
information about the position and energy of the indi- 
vidual particles in the hadron beam, an amplifier which 
takes this signal and amplifies it, and then a kicker which 
takes this information and uses it to deliver an energy- 
dependent non-conservative kick which decreases the lon- 
gitudinal energy spread of the hadron beam. 
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FIG. 1: Schematic of Coherent Electron Cooling 
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For CeC, the pick-up is a co-moving electron bunch 
and hadron beam in a drift, where the individual hadron 
signals are the Debye screened charge perturbations de- 
scribed in Q. The amplifier of the signal is the free- 
electron laser which we describe in this paper. The kicker 
is a chicane which offsets the hadrons from their initial 
signals so that they are displaced from a local maxi- 
mum of the electron density in such a way that hadrons 
with energy greater than the design energy lose energy, 
whereas hadrons with energy less than the design energy 
gain energy. In both the kicker and the pick-up, it is nec- 
essary that the time co-moving between the hadrons and 
the electrons be shorter than a full plasma oscillation or 
else the signal will be greatly diminished. 

Because the bandwidth of the coherent kicks from the 
amplified signal is on the order of the resonant wave- 
length of the FEL, which for most CcC applications is 
on the order of a few hundred nanometers, the cross- 
coherence that arises in stochastic cooling is negligible. 
Thus, the cooling system will continue to reduce the en- 
ergy spread of the hadron beam until another effect is 
encountered. 

To describe the FEL process, we present a theory that 
follows closely the derivation for the FEL instability de- 
rived in [ll| , with the slight modification that the trans- 
verse profile of the electron beam is uniform. We then 
inject this result into the existing results for the kicker 
and pick-up, and determine an exact form for the cool- 
ing decrement. But first we begin with a single-particle 
description of the dynamics. 



III. MAXWELL- VLASOV EQUATIONS 

Consistent with the derivations of the high gain FEL 
in [T2j and [isj and summarized in [TTj . we begin with 
the equations of motion for the single particles in an un- 
dulator subject to the radiation field generated by the 
collective dynamics of the rest of the beam. The hamil- 
tonian equations of motion for small energy deviation and 
high energy (7 3> 1) are given by 



dt 



dz 
1 



1 



1 



edA, 



dt 
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where po = ^/c, Ayj = B-u, /k^ {cos k^z iy — sin. k^z Cx) 
is the undulator vector potential (here we only consider 
helical undulators), A± is the laser field and A^ is the 
longitudinal space charge. The scalar potential has been 
removed by choice of gauge transformation, and Pz has 
been used as the generator of longitudinal translations. 

The Vlasov equation is derived from the conservation 
of single particle phase space volume, so that 

di__dj_ <mdl_ (Udl 
dz dz dz 8% dz dt 







(2) 



Following along with the canonical description of insta- 
bilities in plasmas [l^ . we assume that the phase space 
density of the electron beam is given by / = /i -f- /o, 
where /o is a thermal background and /i is the insta- 
bility. Furthermore, we assume that |/i| ^ |/o|. This 
justifies (i) dropping the term proportional to A\ that 
would appear in equation ()lbp and (ii) dropping terms 
proportional to or higher. Carrying these approxima- 
tions out and knowing that A±_ , oc /i , we obtain the 
equation of motion given by 



df 



Bz ' cil+2^'^ 




(3) 



= 



where K = eAw/rrieC^ is the undulator parameter. Ab- 
sent from this description is an accounting for the trans- 
verse betatron oscillations that arise from the confining 
FODO lattice used on the electron beam in the undula- 
tor. In fact, all the transverse dynamics of this theory 
arise from the Maxwell equations, and it is assumed that 
the current distribution will follow this transverse distri- 
bution. 

By solving the single-particle equations of motion for 
an electron in an undulator, this leads to the relationship 



K ( cosfc,,,2 



70 



3z 



(4) 



where jz ~ — ec / dHfiiH, z,t) is the longitudinal cur- 
rent density. We consider the transverse laser field in 
Fourier space, where its Fourier transform is defined by 



V2^ 



diy d^k^ e'^^-''^e"'^'-^'^^-'^A^{z,i^,kj_) 

(5) 

The transverse Maxwell equation, when Fourier trans- 
formed over ri, is given by 



1 



2tt 



dvd^k I ( — fc^ 



C 



(6) 



It is assumed that the envelope function A±_ is slow- 
varying in the longitudinal direction, and so higher order 
derivatives in z arc small compared to the first deriva- 
tive. This allows us to drop terms that go as over 
kfdz Aj_ . 

By dropping oscillating terms that are 2fc^z out of 
phase with the laser field and defining the Fourier trans- 
form on jz by 



1 



Jz 



/27r 



df d^k±e"'^-'^^e 



(7) 



ikuZ —%k_^cz/2i'ujj. 

^ ^ JZ 
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we obtain for A^, ■ A± in Fourier space the expression 



ittK 



ViOrlO Jo 



%dz' (8) 



where the initial laser field has been set to zero as is the 
case for CeC. 

For the proof of principle, space charge will be a non- 
negligible component of the system. To account for space 
charge, we consider the longitudinal electric field given by 



dtE^ = 

c 

which, under this Fourier transform, gives 



~ Am ~ 

Ey, = 



(9) 



(10) 



All of this is identical to the one-dimensional theory in 
[m except the additional phase factor of k\cz /2vujr that 
appears in the definition of jz, which acts as a detuning. 

By applying an identical Fourier transform of the type 
performed on the current density to the phase space den- 
sity, and assuming that the thermal background is given 

by 



/o = noF{E)R{r^) 



we obtain the coupled Maxwell- Vlasov equation for the 
phase space density of the FEL amplified electron bunch 
with an initial phase space perturbation: 



/l(i,zy, fcx) Q-^(.k^{'^-^)+'^k„y£/eo-kic/2,,u,)z \o + [ dz' e*('=™(l-'')+2/c„>^£/£o-fcic/2i/t^,)(2:'-z)^ 

Jo 

d^q - ) _L / d,J' _ e j^(g-) \ no—R{q- k^) 

toC \ VOJrlO Jo I CVUJr dt 



where F = F{£) is the normalized energy distribution 
and H ~ £ + So where So is the average energy of the 
electron beam. R is the Fourier transform of the trans- 
verse bunch profile. 

The equation of motion is identical in form to that of 
the one-dimensional theory in [Tl| , with the exception of 
the added transverse detuning term k\c/2viUr- Regard- 
less of whether the beam is infinite or finite in transverse 
extent, the inverse gain length, given by 



\2'!ive^ KkwTio 



-1/3 



and Fierce parameter, given by 

P = rfc~^ 

are unaffected by the three-dimensional effects. 

To obtain the longitudinal current density, we take the 
definition jz ~ — ec J d£ fi to equation (jlip . Introducing 
the normalized detuning, space charge parameter, energy 
and transverse wave vector as 



C={l-u)/p 



S-Ke^noV-^{l + K^) 



£ = 2v£/ p£o 



(12a) 



(12b) 



(12c) 



P = k^cT-^/2vujr (12d) 
gives the cleaner and dimensionless form 

Uz,C,k^) ^ -ec^ / d£ e^ic+i-kl),f^ ^ 



(13) 



dF 



dz" Jz{q) + iKliM -jR{q- fc±) 



d£ 



At this point, the method of solution depends on 
whether the beam is to be considered finite or infinite 
in transverse size, which is to say whether the transverse 
dimension of the electron bunch rg is large compared to 
the diffraction length scale of the FEL, d = cT~^ /2vujr- 

On the other hand, whether the transverse spacial ex- 
tent of the initial perturbation can be modeled profitably 
as a delta function in real space (which would be much 
simpler) depends on the comparison of the Debye radius 
to the transverse length scale, ro/d. If the Debye ra- 
dius is much smaller than d, r^/d 1, then the physics 
of a point-perturbation in transverse space should match 
very closely the physics of the initial phase space pertur- 
bation, li r£i/d ^ 1 then the actual physical distribution 
is necessary. If rp/d ^ 1 then we expect the FEL to 
be essentially one-dimensional. By necessity, ro ro 
for the models utilized in Q and @ to be valid. These 
considerations hold for both the infinite and finite beam 
solutions. 
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IV. INFINITE BEAM SIZE 

We first consider a beam that is infinite in tlie trans- 
verse direction, as it is analytically simpler than the fi- 
nite beam size but still contains a reasonable amount of 
physics in its own right. This can be considered in terms 
of the ratio tq / d, where rg is the typical transverse width 
scale of the electron beam and d is the diffraction length 
scale of the FEL. liro/d ^ 1 then the beam is effectively 
infinite and the treatment in this section is useful. Oth- 
erwise the finite beam solution of the next section needs 
to be employed. 

For an infinite beam, R{q — kj_) = 6{q — k±) so the 
above equation (fTTj) reduces to 

I = -ec^ I d£ e^«5+^-fei)v; \,^^ + 




(14) 



This is identical in form to the equations of motion for 
the one-dimcnsional FEL 



ll| with the identification of 



CsD ^ C ^ k^- Due to this similarity, we omit many of 
the details and cut to the solution by Laplace transform 
for the current, which is given by 



s+i(CaD+£) 



fl \z=0 



s- D{1~ isKl) 



where 



D 



d£ 



dF 



1 



d£ s + i(C3D + £) 



(15) 



(16) 



determines the dispersion relation. Equation ()15p gives 
immediately the linear response function in Laplace space 



for the current density perturbation versus an initial 
phase space perturbation 



= / d£'JC{s,C:iD.£')~h lo(f') 



(17) 



such that 



1 



1 



s-D{l- isAD s + liCsD + £') 

(18) 



JC is the linear response function of the modulated cur- 
rent density to an initial phase space perturbation. We 
will use this function to calculate a Green function for 
the FEL phase space distribution, which we will denote 
Gfel- 

By inserting equation back into equation (fTTj) . and 
Laplace transforming for /i in the z coordinate, we ob- 
tain a comparable expression to eqn. ()15|1 for the phase 
space density of the perturbation of the e-beam given by: 



fi{s,C3D,k±,£) 
1 



1 



S + l{C3D + £) 



s + i(C3i5 + £) 
1 



/i lo {£,C:iD.ki_) + 
1 



s- D{1- isAj) ' s- D{1- isAl) 

1 



dF P - ~ - - 
-J / d£'h lo i£',C3D,ks 
d£ J 



S + t{C3D + £') 
(19) 



The form of this equation allows us to write down the 
Green function for the phase space density of an infinitely 
wide e-beam in an FEL amplifier as 



1 



S + liCsD + £) 



Qfel{s, C3D, k±,£;£') 



+ tAts- 



1 



S + t{C3D + £) 



5{£-£') + 



D{1 - isAD 



s- D{1~ isAl) 



dF 



(20) 



d£ S + l{C3D + £') 



where the new FEL phase space density in Laplace- 
Fourier space is given by 

/i(s, CsD, k^,£) = 
f . (21) 

J d£' gFELis,C3D,k^,£\£')MC3D,k^,£')\0 

It is interesting to note that this Green function can be 
clearly divided into two parts. The first part represents 



Landau damping and single-particle non-cooperative mo- 
tion in the FEL undulator. This process does not lead 
to gain, and the term representing it can be dropped in 
a description of the FEL process. The second part con- 
tains the growing roots of the dispersion relation, and 
represents the cooperative gain process of the FEL. It is 
this Green function that is of practical application for the 
theory of Coherent Electron Cooling. 
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The dynamics in the z variable are determined by the 
roots of the dispersion relation, given by 



D 



1 



= 



(22) 



There is another pole from the s-\-i{C^d ~'r£) term in the 
denominator, but the pole associated with this term will 
either oscillate or decay, and therefore does not represent 
amplification as a result of the FEL process, but rather 
a Landau damping of the initial perturbation due to its 
own energy spread. 

As an example calculation, we consider an initial phase 
space perturbation that is monoenergetic, instantaneous 
in time, and a point source. We place this in the context 
of a cold electron beam, where the dispersion relation is 
well known. 

In Fourier space, the transform of the initial condition 
is given by 



(23) 



where it is infinitely broad in the k±_ and C variables. 

Inserting this directly into the Green function calcula- 
tion and taking the inverse Laplace transform on s gives 
a sum with three purely oscillating terms and with the 
three modes of the FEL process. The resulting expres- 
sion is extremely cumbersome, and its physical intuition 
is embodied already in the Green function. We therefore 
only consider the single growing root of the FEL process 
from here on. 

The phase space density is then approximately given 

by 

/l(^,C'3I.) ^ ^ 



S+ + z(C3D +S)s+ + l{C3D + £o) 



1 + iAls+ 



l-D' 



dF 

-J exp (s+z) 
at 



(24) 



\s+ {1 - is+AD - D \ 

where s_|_ is the root of the dispersion relation with 
positive real value. Expectedly, all dependence on fcj^ 
has dropped out, and only Csd remains as the natural 
Fourier parameter for the infinite electron beam. 

Recall the definition of the Fourier transformed phase 
space density as: 



271- 



dv d^k^e'^^-^^: 



-/2-"-/i(z,^.,fc^,£) 



(25) 



It would now be useful to transform the integrals into 
integrals over C^d and k±_ in order that we can determine 
the dynamics of this initial perturbation in real space. 
Recalling the definitions of the parameters leaves 

/i(z,C,r±,f) = ^^^e'« I dC^n d^k 



V2^' cpF- 



(26) 



'fliC3D,Z,£) +C.C. 



where ^ = ujr{z/c — t) +kuZ is the ponderomotive phase. 
It is interesting to note that, although the detailed tem- 
poral information cannot be extracted from this integral 
immediately, the transverse profile can be calculated di- 
rectly as 



/i oc 



1 



2 + K 



-"■l/(4(pC+i) 



(27) 



Because this is a pure phase, it has no transverse size 
information intrinsic to it. The trouble arises from the 
equal value given to all k± by an infinitely small point 
source, which allows the signal to propagate transversely 
instantly as we have not properly accounted for a Lorentz 
covariant description of the transverse electron dynamics. 

A slightly less mathematically pathological case is to 
consider an initially Gaussian transverse distribution, in- 
finitely short. In this case, the previous separation also 
occurs and the resulting width goes as cr? ~ (T§ — i{p£, + z) 
and the profile is gaussian rather than sinusoidal. 



V. FINITE BEAM SIZE 

Having considered the simpler case of the transversely 
infinite beam, we now turn our attention to the case of 
a finite transverse beam profile. To achieve this, we con- 
sider an expansion in the eigenmodes of the transverse 
beam profile, as the Maxwell- Vlasov equation for a finite 
beam is an integral equation with the beam profile func- 
tion as its kernel. From there, we can separate out the 
transverse and longitudinal dynamics, and observe that 
in real space there is no spreading of the eigenmodes, 
consistent with optical guiding. 



A. Eigenmode Expansion 

For the case when R is not a delta function, it is benefi- 
cial to expand the current density solutions in the eigen- 
modes of the R kernel defined by 



1 



d^qR{k-q)M<l) 



(28) 



where for this section we drop the overhats and subscripts 
to simplify the notation. This is best calculated by ex- 
panding R{k — g) as a matrix in terms of some orthonor- 
mal basis. We shall consider such an example calculation 
later, but for now we assume such an eigenbasis is already 
known. 

For any reasonably smooth definition of the transverse 
beam profile, R{k — q)= R{q— k), that is that the kernel 
of the eigenvalue equation is hermitian [l^. This being 
the case, we know that the eigenvectors are orthogonal 
and the eigenvalues are all real, with the orthogonality 
condition being 



d^qd^k ipiik)ipmiq) = Sim 



(29) 
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given that the eigenfunctions are square normalized. 

Expanding the integral of the longitudinal current den- 
sity in a series of the eigenmodes gives 



dz' l{z')^Y.Mk)±e'''''ae{z) 



(30) 



Looking back at the definition of the Fourier transform 
for the current, it is clear that the e'^^^^ terms will cancel, 
and there is no change in the transverse extent of the 
current perturbation, which is consistent with the optical 
guiding discussed in the literature [l6l |. 

The current equation (fH)) can then be reduced to a sys- 
tem of coupled equations for the expansion coefficients. 
That equation is given by 



-ec- 



2v 



dE / d^ki 



d£ 



dz' e*(6+£)(i'-i)x 



(31) 



an + ikl [a'^ + iQ 



dF 



where Qm,e = J d'^k k'^ipm{k)'ipi(k) measures the cou- 
pling between the different eigenmodes of the transverse 
beam profile. This system of equations may be solved 
by Laplace transform (see the Appendix for details) and 
leads to the equation 



(1 + iklu„,)Qi., 



(32) 



where /f is the ifji component of the initial phase space 
perturbation, and 



dE 



dF 



1 



dE s + i{C + E) 



(33) 



determines the dispersion relation for each growing mode. 

It is worth noting, at this point, that the £ index could 
refer to multiple indices, particularly since this is a two- 
dimensional model it could refer to both the azimuthal 
and axial indices, as will be the case when we consider 
the Gaussian beam profile below. For that particular 
case, the different azimuthal modes are uncoupled in the 
Q matrix, so the radial modes for a particular azimuthal 
mode are the ones coupled by the Q matrix, while dif- 
fering azimuthal modes do not mix. This will become 
apparent during the calculation below. 



B. Gaussian Profile 

As an example of this calculation, we consider a Gaus- 
sian transverse beam profile. The procedure for solving 
the initial value problem is as follows: 



1. Calculate the eigenfunctions and corresponding 
eigenvalues to equation psp 



2. Calculate Qm,f. to determine the correct dispersion 
relation 

3. Invert equation (j32p and solve for the initial value 
problem 

Each of these steps should be identical for any other 
transverse bunch profiles; we present only the Gaussian 
case here. 

We begin with a Gaussian beam profile, whose Fourier 
transform is given by 



R{k^ -q) = 



L 



exp ■ 



(fc± - q)' 



and the eigenfunctions therefore satisfy the equation 

'ipi{k±) = 
^ {k±-q)' 



(34) 



- / d'q 



/2^ 



exp ■ 



2L- 



(35) 



It is most convenient to consider this particular form in 
Cartesian coordinates, and in keeping with this we ex- 
pand 



V'f(p) = Xm{Px)Xn{Py) 



(36) 



where each of the individual x satisfy an eigenvalue equa- 
tion of the form 



XmiPt 



1 



dp: 



L 



exp 



{-{p' +P? - 2p,p[)/2L-^] xM 



(37) 



where the resulting eigenvalue for is given by uji = 
A„A,„. It is convenient to define the normalized variable 
/i ~ PiL so that the above eigenvalue equation is given 
by 

1 f°° 

Xm(Ai) = J— + m'^ - 2/i/x')/2} Xm(Ai') 

A„i J-oo 

(38) 

where A„j = A,„V27r. The appropriate scaling for the 
transverse beam size for the full eigenvalue is given by 



2^ 



where uii = Am An. To calculate the normalized eigenval- 
ues, we expand the kernel of this single- variable integral 
equation in terms of Hermite polynomials, as they are 
already related to the paraxial Maxwell equations [17 1. 

It turns out from the properties of Hermite polynomi- 
als that only the evens and odds couple, so each Xm is a 
series in either even or odd Hermite polynomials. In this 



7 



case, the matrix equation for the even Hermite polyno- 
mials is given approximately by the matrix elements 

dfi / d^i' exp {-(^i^ + - 2/i^')/2} 

-oo J — oo 

(39) 

Furthermore, to good approximation, the expansion 
can be carried out for the first two Hermite functions in 
the series. We therefore consider the two-mode case. For 
the principle even mode, the matrix is given by 



G = 



9 1. / ZZL 

^ V 3 3 V 3 



(40) 



for the vector components (7Jo(m)7 ^2(m))* exp(— ^^/2). 
The eigensystem here has eigenvalue Xeven = 2.2382 with 
corresponding eigenvector 



and a smaller eigenvalue A2 
eigenvector 



.9294 
.3690 



.83178 with corresponding 



-.1465 
.3690 



To validate these numerical results we take the matrix 
to next order, i.e. to order Hi{^) in the expansion, and 
the matrix is given by 



2a/¥ 



17 
54^ 



2-jT \_ \ 



1 

9 

17 
54 
227 
324 



(41) 



which yields an eigensystem given by X'^ = 2.3157, A2 = 
1.2005 and A3 = .27073 with corresponding normalized 
eigenvectors 



Vl 




V2 




matrix expansion, and from analysis of the eigenvector 
components the H^^fi) level of expansion is negligibly 
small compared to the other two components for the 
eigenvector with the maximal eigenvalue. 

Carrying out a similar procedure for the Hi (/i) - H3 (/i) 
eigenmode gives a maximal eigenvalue Xodd = 1.7161 and 
eigenvector 



Vodd 



.8456 
.5339 



It is now necessary to calculate the various matrix ele- 
ments for Q. For the purposes orderly book-keeping, we 
define the following modes 



)Ceven everi (My) 



Ipodd — Xodd{lJ-x)Xodd{^^y) 



(42a) 



(42b) 



V2 



iXoddifJ-x)XevenilJ-y) + Xe 



{lJ.x)Xodd{fiy)) 

(42c) 



{Xodd{^.x)Xe 



i(A*y) ~ Xeven{lJ-x)Xodd{iJ-y)) 



(42d) 

as the orthonormal basis of expansion. The correspond- 
ing eigenvalues are given by uJeven ~ 5.0095, Cjodd = 2.945 
and ujj^ = w_ = 3.8410. Under this particular basis the 
Hermite polynomials have a particularly nice relation for 
the Q matrix elements, and Q is diagonal. The indi- 
vidual modes do not couple, and their growth rates are 
determined by the dispersion relation 



Duj„ 



(1 -f isA;)] + (1 + iA; )Q 



(43) 



The individual Q are given by Qeven ~ 2.5144:6/ L^, 
Qodd = 6.35275/i'^, and Q+ = g_ = 4.43333/L'^. The 
growth rate for these parameters is given in figure ([2|). 
with L = 3. 

To recap, we have calculated an eigenbasis for the 
transverse beam profile, yielding a linear superposition of 
even- and odd-numbered Hermite polynomials, and their 
corresponding eigenvalues. The series is truncated at two 
dominant modes, and because of the particular nature of 
the Hermite polynomial expansion basis, the Q matrix 
is diagonal. If Q had off-diagonal matrix elements, there 
would be "gain leakage" between the connected eigenvec- 
tors. 




.03343 \ 
-.1879 
.2245 J 



We can conclude from this that the largest eigenvalue 
can be accurately determined to within 3% with the 2x2 



C. One-Dimensional Limit 

Because the eigenvalues are totally independent of the 
transverse size, and only Q is dependent, it is straightfor- 
ward to get directly to the one-dimensional beam limit 
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FIG. 2: Growth rates for three eigenmodes: (i) top is of mode 
with largest eigenvalue, (ii) is degenerate case of the odd/even 
mixtures, (iii) is of smallest eigenvalue 



for the dispersion relation. By redefining the normaliza- 
tion as 

S = sLo-y^ (44a) 
C - Cc.^i/s (44b) 

= kyj' (44c) 

Qra = OmC^™'/' (44d) 

the dispersion relation takes the form 

' - i~ %n2 ^^ + + + ^K^T)Q = (45) 

The actual scaling is such that, for large beams, the 
portion of this dispersion relation identical in form to 
the one-dimensional dispersion relation comes to strongly 
dominate over the perturbation correction for finite size, 
taken by the value of Qm- For the case of an infinitely 
large transverse size all functions are eigenmodes and all 
all eigenvalues are unity, therefore we can obtain the one- 
dimensional limit through this limit. 

VI. DISCUSSION 

We have presented a theoretical model for the dy- 
namics of a high-gain free-electron laser with three- 
dimensional effects. The model is analytically solvable 
up to a numerical Fourier transform, and for that reason 
is useful for benchmarking the massive tracking programs 



used to simulate FELs. All results in this paper are re- 
duced to a handful of dimensionless numerical Fourier 
transforms. 

When applying the finite beam case, we observe that 
only the principle four modes grow rapidly. The higher 
order modes have eigenvalues substantially smaller than 
these modes, and can be neglected in comparison to the 
principles. We can therefore conclude from this model 
that an FEL can be effectively characterized by only a 
handful of well-understood eigenmodes. Furthermore, 
this particular model includes optical guiding by con- 
sideration of the transverse eigenmodes of a stationary 
beam. By contrast, we observe spreading of the infinite 
beam case at a slower than linear rate. 

The principle goal of this solution to the three- 
dimensional FEL equations is to develop an understand- 
ing of the charge modulation at the end of the undula- 
tor. A thorough understanding of the phase information 
of the FEL instability is necessary to properly calibrate 
the chicane and inject the hadrons with a proper dis- 
placement with respect to the local charge maxima of 
the bunch. This model provides the phase information 
up to a three-dimensional Fourier integral, which is well- 
bounded and provides adequate benchmarking for nu- 
merical simulations. 

The existing analytical models for the kicker and pick- 
up of CeC involve an infinitely large electron beam, or 
equivalently that the initial perturbation be small com- 
pared to the transverse size of the electron beam. The 
results are also obtained analytically for the k — 2 dis- 
tribution. To match up with these theories, we consider 
the case where R(k±^ ^ <?) = 5{k±^ — q) and with the corre- 
sponding dispersion relation for a k — 2 distribution. The 
results for k — 2 are not presented in this paper, but it is 
straightforward to obtain the dispersion relation from the 
dispersion integral, and we can now consider a complete 
description of the phase space evolution of the electron 
bunch through the CeC process. 

This analytical model was developed to provide bench- 
marking for the proof of principle CeC system to be 
implemented at RHIC. For the FEL for the proof of 
principle, the transverse size of the electron bunch is 
To ~ 3 mm, the resonant wavelength is ~ .5 /im, 
and a gain length of approximately F^^ = 3 m. In this 
case, the transverse length scale d « .35 mm and it is ex- 
pected that the three-dimensional infinite beam theory 
should be a reasonable description of the FEL amplifier 
portion of CeC. 

At present this model has no way of coping with a 
transverse momentum spread in the initial phase space 
perturbation or with betatron oscillations, because all of 
the dynamics are taken directly from Maxwell's equa- 
tions. As such, it is not clear what effect transverse mo- 
mentum spread and betatron oscillations will have on 
the phase information of the amplified signal. Numerical 
modeling or a more complete theoretical description are 
necessary to account for these effects. 
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